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Abstract. The production at central rapidity of Kg, A, 5 and Vl particles in Pb-Pb 
collisions at 158 A GeV/c has been measured by the NA57 experiment over a centrality 
range corresponding to the most central 53% of the inelastic Pb-Pb cross section. In 
this paper we present the rapidity distribution of each particle in the central rapidity 
unit as a function of the event centrality. The distributions are analyzed based on 
hydrodynamical models of the collisions. 
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1. Introduction 

Lattice quantum chromodynamic calculations predict a new state of matter of 
deconfined quark and gluons (quark gluon plasma, QGP) at an energy density exceeding 
~ 1 GeV/fm 3 1J. Nuclear matter at high energy density has been extensively 
studied through ultra-relativistic heavy ion collisions (for recent developements, see 
reference [2]). 

Within the experimental programme with heavy-ion beams at CERN SPS, NA57 
is a dedicated experiment for the study of the production of strange and multi-strange 
particles in Pb-Pb collisions at mid-rapidity 

The measurement of strange particle production provides one of the most powerful 
tools to study the dynamics of the reaction. In particular, an enhanced production of 
strange particles in nucleus-nucleus collisions with respect to proton-induced reactions 
was suggested long ago as a possible signature of the phase transition from colour 
confined hadronic matter to a QGP [3]. The enhancement is expected to increase with 
the strangeness content of the hyperon. These features were first observed by the WA97 
experiment [5] and subsequently confirmed and studied in more detail by the NA57 
experiment Other insights into the reaction dynamics have been inferred from the 
p T distributions of the strange particles: the study of the transverse expansion of the 
collision and the p T dependence of the nuclear modification factors have been presented, 
respectively, in reference [2j and [B]. 

Rapidity distributions provide a tool to study the longitudinal dynamics; e.g., 
differences between protons and anti-protons have been interpreted as a consequence 
of the nuclear stopping j^j. If hyperons, like protons, keep a 'memory' of the initial 
baryon density, then the relative pattern for the rapidity distribution of hyperons and 
anti- hyperons should resemble that of protons and anti-protons [TU] . 

Hydrodynamical properties of the expanding matter created in heavy ion reactions 
have been discussed by Landau [11] and Bjorken [T2j in theoretical pictures using 
different initial conditions. In both scenarios, thermal equilibrium is quickly achieved 
and the subsequent isentropic expansion is governed by hydrodynamics. 
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2. Data sample 

The NA57 experiment has been described in detail elsewhere Ej- The results 
presented in this paper are based on the analysis of the Pb-Pb data sample collected 
at 158 A GeV/c (460 M events); the sample analyzed here corresponds to the most 
central 53% of the inelastic Pb-Pb cross-section. The data have been divided into five 
centrality classes, labelled as 0,1,2,3 and 4 from the most peripheral to the most central, 
according to the value of the charged particle multiplicity sampled at central rapidity 
by a dedicated silicon microstrip detector. 

The procedure to measure the multiplicity distribution and to determine the 
centrality of the collision for each class is described in reference jT^j. The fractions of 
the inelastic cross-section for each of the five classes, determined assuming an inelastic 
Pb-Pb cross-section of 7.26 barn (as calculated with the Glauber model), are given in 
table Q] 



Table 1. Centrality ranges for the five classes. 



Class 4 


3 


2 1 


ffMneZ (%) tO 4.5 


4.5 to 11 


11 to 23 23 to 40 40 to 53 



3. Analysis and results 

A detailed description of the particle selection procedure, as well as of the corrections for 
geometrical acceptance and for detector and reconstruction inefficiencies, can be found 
in reference |7j. Here we describe the principal features of the analysis. 

The selection of strange particles was accomplished by reconstructing their weak 
decays into final states containing only charged particles: 

Kg — ► 7T+ + 7T~ A — > p + 7T~ 

S" -> A + 7T- Q- -> A + K- ^ 

U P + 71 p + 7X 

and their corresponding charge conjugates for antihyperons. Geometric and kinematic 
constraints were used to obtain high purity strange particle samples: the total amount 
of combinatorial background is estimated to be 0.7%, 0.3% and 1.2% for Kg, A and A 
respectively, and less than 5% for 5 and fl. 

An event by event procedure, using a GEANT-based ^Hj Monte Carlo simulation 
of the apparatus, was applied to correct the data for geometrical acceptance and 
for detector and reconstruction inefficencies. This procedure, which is very accurate 
but rather time consuming, assigns a weight to each reconstructed particle and was 
specifically developed for rare signals such as Q and S. The same procedure was also 
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applied to small fractions of the much more abundant Kg and A samples. The selected 
particles have been sampled uniformly over the whole data taking period; the sizes of 
those subsamples were chosen in order to reach a statistical accuracy better than the 
limits imposed by the systematic errors. 

Acceptance regions have been refined off-line in such a way to excluded from the 
final sample the border regions where the correction factors would have larger systematic 
errors. The selected regions of the transverse momentum vs. rapidity plane are shown 
in figure HJ 
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Figure 1. The y-pi acceptance windows. Dashed lines show the position of mid- 
rapidity (y cm = 2.92). 



Extensive checks were performed by comparing real and Monte Carlo distributions 
for several geometrical and kinematical parameters, in particular those used in the 
particle selection procedure. In all cases a good agreement between MC and data has 
been found [3 HH] . 

The double-differential distribution for each particle species has been parameterized 
using the expression 



d 2 N dN 



p T dp T dy dy 



oc — — ■ exp 



a/ m 2 + 



Pt 



T, 



(2) 



app 
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where the inverse slope parameter T app ('apparent temperature') has been extracted by 
means of a maximum likelihood fit of equationEJto the data jTj. By using equationEJwe 
can extrapolate the rapidity distribution measured in the selected acceptance window 
down to p T = 0: 



The contribution of the extrapolation procedure to the systematic errors on has been 
evaluated using an alternative parameterization for the p T dependence of the invariant 
double-differential distribution: instead of the single exponential in equationEl the blast- 
wave model, a hydrodynamics-inspired model with a kinetic freeze-out temperature Tf 
and a linear transverse flow velocity field f3±(r) [T7|, was used with freeze-out parameters 
Tf = 144 MeV and < (3± >= 0.38 which we determined from the blast- wave analysis 
of the p T spectra [7j. This alternative procedure does not affect the shape of the ^ 
distributions but only the absolute values of the extracted yields; the yield of the Kg 
meson increases by 6%, that of the A hyperon remains unchanged, those of the H and 
Q hyperons decrease by 5% and 20%, respectively. This results in the main source of 
systematic uncertainty for the Q hyperon. 

3.1. Results in the centrality range 0-53% 

The measured rapidity distributions for the centrality range corresponding to the most 
central 53% of the Pb-Pb inelastic cross-section (total sample) are shown in figureEJwith 
closed symbols. For all hyperons the rapidity distributions are found to be symmetric 
with respect to the rapidity of the centre of mass ('mid-rapidity') within the statistical 
errors. A similar conclusion cannot be drawn for Kg since our acceptance coverage does 
not extend to backward rapidity. The symmetry of the Pb-Pb colliding system allows 
us to reflect the rapidity distributions around mid-rapidity. In figure El open symbols 
are used to plot the measured points after such a reflection. 

The rapidity distributions of A, S~, H + and Q are compatible, within the error bars, 
with being flat within the NA57 acceptance window. For the Kg and A spectra, on the 
other hand, we observe a rapidity dependence. The rapidity distributions for these 
particles are well described by Gaussians centred at mid-rapidity. 

The total yields integrated over one unit of rapidity as reported in references [HJ E| 
are also shown in figure El as shaded bands superimposed to the corresponding rapidity 
distributions. In these bands the dotted and full lines represent, respectively, the 
statistical and the systematic errors. These yields are calculated with a maximum 
likelihood fit using the Gaussian parameterization of dN/dy for K s and A, and a flat 
dN/dy for all other particles. The numerical values of the yields are reported in table El 
along with their errors §. 

§ In previous conference proceedings and in reference ^H] the gj-^jp distributions of Kg and A were 
evaluated assuming flat rapidity distributions: the resulting yields are equal to 11.7±0.3 and 1.17±0.03, 
respectively; the difference is less than the statistical error. 




(3) 
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Figure 2. Rapidity distributions of strange particles in the most central 53% of Pb-Pb 
interactions at 158 A GeV/c. Closed symbols are measured data, open symbols are 
measured points reflected around mid-rapidity. The A and 5 results have been scaled 
by factors 4 and 2, respectively, for dispaly purposes. The superimposed boxes show 
the yields measured in one unit of rapidity (as published in references [fil llSj ') with the 
dashed and full lines indicating the statistical and systematic errors, respectively. 



Table 2. Absolute yields in one unit of rapidity calculated by means of maximum 
likelihood fits. The first error is statistical, the second is systematic. Values of the 
X 2 /ndf for the measured dN /dy distributions calculated assuming dN /dy constant and 
equal to the absolute yields are also quoted. 





JUcm—0.5 dy & 


X 2 /ndf 




11.8 ±0.3 ±1.2 


29/6 


A 


1.14 ±0.03 ±0.12 


12/6 


A 


7.84 ±0.21 ±0.78 


3.9/6 




0.852 ±0.015 ±0.085 


3.7/7 




0.209 ±0.007 ±0.021 


3.6/5 




0.172 ±0.013 ±0.026 


5.0/4 



The Gaussian parameters which minimize the x 2 of the Kg and A distributions are 
given in table El 
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Table 3. Results of fits with a Gaussian to the rapidity distributions of Kg and A. 





dN 
dy 


y=yc m ° X 2 /ndf 


K° 
A 


13.2 ±0.4 0.69 ±0.07 2.5/5 
1.22 ±0.04 0.83 ±0.22 3.0/5 



3.2. Comparison with results 

The large acceptance NA49 experiment has measured the rapidity distributions of a 
number of particle species over a broad rapidity range (« 3 units of rapidity centred 
around mid-rapidity) in central Pb-Pb collisions ^H]- The A rapidity distribution 
measured by NA49 is flat over the experiment's ^/-coverage and the A distribution has 
been fitted to a Gaussian with a width of 0.95 ±0.05 [20J, compatible with our findings. 
NA49 has fitted the rapidity distributions of negatively and positively charged kaons 
with a two-Gaussian parameterization [21]. Assuming the Kg to be an average of K~ 
and K + , and taking the (K~+K + )/2 average of the NA49 fitted parameters one obtains 
in our rapidity range (2.93 < y < 3.70) rms = 0.219 to be compared with rms = 0.221 
of our measured distribution. For the multi-strange 5 and Q hyperons, our rapidity 
distributions are compatible within the error bars both with being flat and with the 
shapes of the NA49 fits [221 123] • 

In conclusion, there is a fairly good agreement between the two experiments as far 
as the shape of the rapidity distributions are concerned. On the other hand a systematic 
discrepancy on the values of the yields is found relative to the NA49 results, which are 
lower for all particles by about 25-30% [21]. The origin of such a discrepancy has not 
yet been understood despite the joint efforts of the two collaborations. 

3.3. Centrality dependence 

In figures [3] and [3] we show the centrality dependence of the rapidity distributions for 
Kg and A, respectively. As in figure [21 the shaded bands represent our determination of 
the yield in one unit of rapidity, according to the maximum likelihood fit [UJEIj with 
the dotted and full lines indicating the statistical and systematic errors, respectively. 
Gaussian y 2 fits to the spectra are also superimposed, with the fit parameters reported 
in table 0] For both particles, the width of the rapidity distributions is constant within 
the errors in the five centrality classes (i.e. from 40-53% to 0-4.5%, see table E}. 

The rapidity distributions in the five centrality classes of A, S~ and S + are shown 
in figures M and respectively. For statistical reasons the rapidity distributions of 
f2~+f2 + have been calculated in three centrality classes instead of five: + 1, 2 and 3 + 4 
and they are shown in figure [S] In all the centrality classes, the rapidity distribution 
of the A hyperon is consistent with being flat over the considered range. In the same 
rapidity range, the A distribution varies by about 40% (class 4). It is likely that the A 
hyperon rapidity distribution reflects the overall net baryon number distribution. The 
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Figure 3. Rapidity distributions of Kg for the five centrality classes. 
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Figure 4. Rapidity distributions of A for the five centrality classes. 
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Figure 5. Rapidity distributions of A for the five centrality classes. 
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same behaviour was observed for the y distribution of protons in central Pb-Pb collisions 
at the same energy by the NA49 experiment "H"]. The 5~ distributions are found to be 
flat in one unit of rapidity like the A distributions. For the rare 5 and f2 particles, the 
NA57 collected statistics do not allow us to observe significant deviations from a flat 
distribution in our limited rapidity range. 



Table 4. Results of fits with a Gaussian to the rapidity distributions of Kg and A for 
the five centrality classes defined in table"]] 





4 


3 
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1 







dN 1 

dy 1 y = Hcm, 


30.3 ±2.2 


24.4 ± 1.6 


17.2 ±0.9 


9.4 ±0.6 


3.3 ±0.4 




a 


0.61 ±0.12 


0.67 ±0.13 


0.72 ±0.13 


0.61 ±0.11 


0.71 ±0.34 




X 2 /ndf 


1.2/4 


2.4/3 


5.2/4 


1.9/5 


2.4/4 


A 


dN 1 

dy 1 y = Ucm 


2.72 ±0.20 


1.90 ±0.13 


1.68 ±0.09 


0.90 ±0.06 


0.44 ±0.05 




a 


0.52 ±0.14 


0.73 ±0.34 


0.9 ±0.5 


0.61 ±0.20 


0.48 ±0.19 




X 2 /ndf 


2.6/4 


2.4/4 


1.14/3 


1.18/4 


7.0/5 



The values of the yield in one unit of rapidity around mid-rapidity as a function of 
the centrality of the collision are given in table [5] for all the particles. 



Table 5. Strange particle yields in one unit of rapidity around mid-rapidity, 
Iy Cm -0 5 7^7 dy, for the five centrality classes defined in table [J Statistical (first) 
and systematic (second) errors are also quoted. 





4 3 


2 


1 


K° 


26.0 ±1.7 ±2.6 21.6 ±1.2 ±2.2 


15.3 ± 0.8 ± 1.6 


8.1 ±0.5 ±0.8 3.05 ±0.33 ±0.31 


A 


2.44 ±0.14 ±2.4 1.81 ±0.10 ±0.18 


1.58 ±0.07 ±0.16 


0.80 ±0.04 ±0.08 0.41 ±0.03 ±0.04 


A 


18.5 ± 1.1 ± 1.8 15.0 ±0.8 ±1.5 


9.5 ±0.5 ±0.9 


5.19 ±0.29 ±0.51 2.30 ±0.22 ±0.23 




2.08 ±0.09 ±0.21 1.80 ±0.07 ±0.18 


1.07 ±0.04 ±0.11 


0.52 ±0.02 ± 0.05 0.181 ±0.013 ±0.018 




0.51 ±0.04 ± 0.05 0.37 ±0.03 ±0.04 


0.29 ±0.02 ±0.03 


0.14 ±0.01 ±0.01 0.045 ±0.007 ±0.004 


Q 


0.39 ±0.04 ±0.06 


0.23 ±0.03 ±0.03 


0.069 ±0.012 ±0.010 



4. Longitudinal dynamics 

The transverse dynamics of the collisions have been studied in reference [7j from 
the analysis of the transverse momentum distributions of the strange particles in the 
frame- work of the blast- wave model [I7|. The rapidity distributions can be used to 
extract information about the longitudinal dynamics. We use an approach outlined 
in reference [T7] (i.e., the same blast- wave model used for the study of the transverse 
dynamics) and [23], where, respectively, Bjorken and Landau hydrodynamics [1*2*1 
are folded with a thermal distribution of the fluid elements. 

In figure [HI the observed rapidity distributions are compared with the expectation 
for a stationary thermal source and with a longitudinally boost-invariant superposition 
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Figure 6. Rapidity distributions of 5 for the five centrality classes. 
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Figure 7. Rapidity distributions of S for the five centrality classes. 
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Figure 8. Rapidity distributions of £1 +f2 for central (0-11%), semi-central (11- 
23%) and semi-peripheral (23-53%) collisions. 
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of multiple isotropic, locally-thermalized sources (i.e. Bjorken hydrodynamics). 
Each locally thermalized source is modelled by an m T -integrated Maxwell-Boltzmann 
distribution with the rapidity dependence of the energy, E = m T cosh(?7) explicitly 
included 

l2 ™ 2 2 , . m . 

rp2 ' T ~1T + T2" ex P -T^coshr] , 4 

Tj Tf cosh 7] cosh rj T f 



where 7/ is the freeze-out temperature and m T = vPt + m2 • The distributions are 
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Figure 9. Rapidity distributions of strange particles for the centrality range 
corresponding to the most central 53% of the inelastic Pb-Pb cross-section as compared 
to the thermal model calculation of equation 0] (dotted lines, in red) and a thermal 
model with longitudinal flow (full lines, in black). 



integrated over source element rapidity to extract the maximum longitudinal flow, rj, 
dN _ f Vmax dN t 
dy 



nit 



Ijrna x 



drj 



(y - v) dv 



Pl = tanh?7 ri 



(5) 



where j3i is the maximum longitudinal velocity in units of c. The average longitudinal 
flow velocity is evaluated as < /3l >= tanh(^ max /2). A simultaneous fit of the 
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function defined by equation to the rapidity distributions of the strange particles 
gives < (3l >= 0.42 ± 0.03 with x 2 jndf = 28.2/32. The freeze-out temperature has 
been fixed to the value Tf = 144 MeV obtained, for the most central 53% of the 
inelastic Pb-Pb collisions, from the analysis of the transverse mass spectra of the same 
group of particles [Zj. In the same analysis the average transverse flow velocity has been 
determined to be < j3± >= 0.38 ± 0.02, i.e. slightly less than the longitudinal velocity 
determined in this analysis; this indicates substantial stopping of the incoming nuclei 
as a consequence of the collision. 

In principle, also the freeze-out temperature can be fitted from the rapidity 
distributions along with the longitudinal velocity. However, the sensitivity to the freeze- 
out temperature is very limited. For instance, changing Tf from 144 to 120 MeV results 
in only a 2% increase of < /3l >■ The partial contributions to the total x 2 from the 
individual particle rapidity spectra are given in table El Within our uncertanties, we do 
not observe any particle to deviate from the common description given by a collective 
longitudinal flow superimposed to the thermal motion. A combined fit performed 
only to the Kg and A rapidity distributions yields a smaller value of the flow, i.e. 
< (3 L >~ 0.36 ± 0.03. It is worth noting that the flattening of the rapidity spectra 
with increasing particle mass, which is also observed in the data, is due in the model 
to the collective dynamics: all particles are driven by the flow with the same velocity 
independently of their masses. 

Table 6. Partial contributions to the \ 2 01 the best-fit to the thermal model with 
longitudinal flow, equation^ The x 2 /ndf of the global fit is 28.2/33, since there are 
six normalization constants (the factors A in equation |5j) and the parameter 



particle 


n. of points 


x 2 


K° 


7 


5.1 


A 


7 


3.6 


A 


7 


6.7 




8 


4.5 




6 


3.1 


n 


5 


5.2 


tot 


40 


28.2 



In Landau hydrodynamics, the amount of entropy {dS) contained within a (fluid) 
rapidity dr] is given by 11J 



dS 
drj 



7T 



R 2 ls f3c s exp[f3ujf} 



h{q) 



(6) 



where q 



f 



c 2 s V 2 , 



Uf = ln(T//T0), c s is the velocity of sound in the medium, 



T is the initial temperature, rj is the rapidity, R is the radius of the nuclei, 21 is the 
initial length, sq is the initial entropy density, 2(3 = (1 — c 2 )/c 2 and Jo, h are Bessel 
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functions. The quantity ttR 2 Isq is used to normalize the spectra at mid-rapidity. The 
particle rapidity distribution is obtained, as for the Bjorken superpositon of 

the multiplicity density in rapidity space (dN/ dt] oc dS/ drj) with a thermal distribution 
of the fluid elements, 



In the Landau model the width of the rapidity distribution is sensitive to the velocity 
of sound and to the ratio of the freeze-out temparature to the initial temperature. 
While integrating over 77 in equation El the range of rj is treated as a parameter in 
case of Bjorken hydrodynamics; moreover in the Bjorken case (equation |3J) the factor 
dN/dr), which appears in equation has been included in the overall normalization 
factor A since the entropy density dS/drj is independent of the rapidity, in accordance 
with the assumption of boost invariance along the longitudinal direction in the 

case of Landau hydrodynamics, the integration limit is fixed by r] max = —rj min = 
\n(T /T f)/c s || and the multiplicity density in rj space (dN/dr/) is written explicitly 
in the r] integration (equation Ej). Landau hydrodynamics can also reproduce 
simultaneously the distributions for all the strange particles considered (x 2 /ndf ~ 
28/32) but we are not able to put stringent constraints on both the velocity of sound 
and the ratio Tj/Tq. The confidence level contours in the c 2 s vs. ^ parameter space 
are shown in figure ITU1 For instance, the hypothesys of a perfect gas (i.e. c 2 s = 1/3) 
would result (at the 3a confidence level) in either Tf/T « or Tf/T ~ 0.6. In fact, 
two physical regions are constrained at the 3a confidence level, the first located at small 
values of Tf/T and the second between 0.5 and 0.8; on the other hand, the region at 
> 1/3 is unphysical. Both physical regions span over the full range < (? s < 1/3. 

5. Conclusions 

We have measured the dN/dy distributions of high purity samples of Kg, A, H and Q 
particles produced at central rapidity in Pb-Pb collisions at 158 A GeV/c over a wide 
centrality range of collision (i.e. the most central 53% of the Pb-Pb inelastic cross- 
section). 

In the unit of rapidity around mid-rapidity covered by NA57, we have performed 
fits to the dN/dy distributions of Kg and A using a Gaussian parameterization: the 
resulting widths are compatible with each other and constant as a function of centrality. 

Contrary to A, the A spectra are flat to good accurancy in the range of rapidity 
and centrality considered; this would indicate that the A hyperon conserves "memory" 
of the initial baryon density. 

The rapidity distributions of the Q particle are found to be flat within the errors 
in one unit of rapidity for central (0-11%) and peripheral (23-53%) collisions. 

|| In reference |26| a modification has been developed (Srivastava) where the integration limit for 
rapidity is infinite, but this case has not been considered in the present analysis. 
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Figure 10. The square of the velocity of sound in the medium (in unit of c 2 ) versus 
the ratio of the freeze-out temperature to the inital temperature. The la (full curves) 
and the 3cr (dashed curves) confidence contours are shown. The dotted line at c 2 = 1/3 
show the ideal gas limit. 



Boost-invariant Bjorken hydrodynamics can describe simultaneously the rapidity 
spectra of all the strange particles under study with x 2 / n df ~ 1, yielding an average 
longitudinal flow velocity < (3l >= 0.42 ± 0.03, sligthly larger than the measured 
transverse flow. The isotropic collective expansion of the system suggests large nuclear 
stopping. 

A fairly good description is also provided by Landau hydrodynamics, which allows 
us to put constraints in the parameter space of the velocity of sound in the medium and 
the ratio of the freeze-out temperature to the initial temperature. 
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